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ABSTRACT 

The analysis of a fast- response fuel- control valve for use in turbojet engine re- 
search is presented. This analysis includes cases where the valve is (1) closely coupled 
to the engine spray nozzles and (2) located some distance from the engine, thus require- 
ing connection by a hydraulic transmission line. The analysis results in a complete non- 
linear model and an equivalent linearized model for both cases. The dynamic perfor- 
mance of the analytical models is presented and compared against corresponding experi- 
mental fuel-flow response data. A discussion is included which demonstrates how oper- 
ating conditions affect fuel-flow frequency response. 
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DYNAMIC PERFORMANCE ANALYSIS OF A FUEL-CONTROL VALVE 
FOR USE IN AIRBREATHING ENGINE RESEARCH 
by Peter G. Batterton and John R. Zeller 
Lewis Research Center 

SUMMARY 

An analysis of a fast-response fuel-control valve with a reducing pressure regulator 
is presented. This type of fuel-control valve is often used for dynamic and steady-state 
engine studies. A nonlinear analytical model for this valve is developed, which accur- 
ately predicts the dynamic performance of fuel flow from the valve in response to an in- 
put command. Since the fuel-flow response can be uncoupled from the valve actuating 
servo, this servo is not considered in this report. Many experimental configurations 
require that the fuel control be located some distance from the engine, thus requiring an 
interconnecting transmission line. The analysis is extended to include a transmission 
line between the valve and the load. 

A simplified linearized equivalent of the nonlinear analytical model is also derived 
and found to be an adequate representation for the valve. This linearized equivalent is 
particularly useful for the prediction of fuel-flow dynamic performance, both bandwidth 
and damping, for any load configuration. Also, the resulting transfer functions are use- 
ful in system performance evaluations. Results are included to demonstrate how this 
linear analytical model can be used to show some of the effects on the dynamic perform- 
ance due to various changes in operating parameters. 

The transmission line degrades the amplitude frequency response of the fuel system. 
In addition, the time delay introduces an undesirable phase shift. Both effects can be 
tolerated for dynamic studies, but they seriously handicap the use of the valve in ad- 
vanced engine control systems. The use of a transmission line should be avoided when- 
ever possible. 


INTRODUCTION 


Renewed interest in turbine engine research studies, both steady-state and dynamic, 



has created a need for a fast-response, high-performance fuel-control valve. Such a de- 
vice should be able to (1) accurately control steady-state engine fuel flow and (2) provide 
controlled dynamic fuel-flow disturbances for the determination of the dynamic character- 
istics of various engine components. With these capabilities, the fuel-control valve can 
be employed within a control loop to facilitate the evaluation of new and more sophisti- 
cated engine control concepts. 

In reference 1, it was concluded that a fuel-control device based on the principle of 
reducing fuel supply pressure to maintain a constant regulated pressure drop across a 
variable area orifice best satisfies these requirements. The hydromechanical portion of 
this device controls output fuel flow in proportion to the area of a fuel metering orifice 
and independent of output pressure. Metering orifice area is controlled by an electro- 
hydraulic servo actuation system. Fuel-control valve performance can be considered 
independently of the performance of the fuel valve servo. Thus, this report is re- 
stricted to fuel-flow performance with respect to a metering orifice area input. 

In order to efficiently use this fuel-control valve for a specific application, it is 
necessary to understand and predict its dynamic performance for each set of operating 
conditions. To facilitate this understanding, we will present a detailed analysis of this 
hydromechanical control device operating directly into an orifice output termination (en- 
gine spray nozzles). Actual experimental installations often require the fuel-control de- 
vice to be located some distance from the spray nozzle load. Therefore, a long line is 
required to connect the two elements, and the dynamics of this transmission line will be 
included in the analysis of the fuel-flow high-frequency dynamic performance. This in- 
clusion will supplement the analysis presented in reference 1. To assist in predictions 
of dynamic performance, we will also derive a simplified linearized analytical model of 
this fuel-control valve and its connected load. 


PHYSICAL DESCRIPTION 

Figure 1 shows the general arrangement of a fuel-control valve with a reducing valve 
pressure regulator. This regulator employs a piston-actuated variable orifice and a pre- 
loaded spring that combine to regulate the pressure drop across the control or fuel me- 
tering orifice. This pressure drop is held to a value proportional to the force of the pre- 
loaded spring. The regulating device is designed to be sensitive enough to maintain this 
pressure nearly constant over the complete range of flows and output pressure variations. 
As a result, the control orifice flow is proportional to the control orifice area. 


2 



High-pressure fuel 
supply input — -** 


c 


Reducing valve orifice — ' 

Reducing valve piston—- 
Reducing valve spring — ■ 


Mv 7 - 


H- 


•R3 


Control 

orifice 

actuator 


■M 


-Control orifice 


/ contre 


, Controlled output 
J fuel flow 


Figure 1. - Fuel-control valve with reducing valve pressure regulator. 


ANALYTICAL DESCRIPTION 

Figure 2 shows a more detailed version of the hydromechanical portions of the fuel 
throttling device of figure 1. It also shows the valve output loaded by a fixed orifice with 
and without the transmission line. This type of load typically represents the spray noz- 
zles of a turbojet engine. In the next sections, a detailed nonlinear analytic representa- 
tion of this hydromechanical device is derived. 
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Orifice Equations 


The basic orifice equation (ref. 2) is 


q 



VF 


ressure drop 


Thus, for the load orifice, equation (1) becomes 


= 




( 1 ) 


( 2 ) 


where 



(3) 


(All symbols are defined in appendix A. ) The pressure p d is the downstream load ori- 
fice pressure. It was held constant at atmospheric pressure for the experimental and 
analytical results presented in this report. In the case of an active engine, it is equal to 
combustion chamber pressure. Its effect might well be included in any system analysis. 
It is included in the equations and computer diagrams of this report for the sake of gen- 
erality. As explained in reference 1, the type of fuel-control valve analyzed in this re- 
port displays the desirable feature of being insensitive to low-frequency downstream 
pressure fluctuations. 

The two variable orifices (control and reducing valve) are rectangular in shape and 
of constant width. Thus, the equations for the flow are as follows: 


where 


^c = K c x c V^r - Pc 



and 


K r x r 




(4) 


(5) 


( 6 ) 
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where 


K = Ch 
r r 


i 


(7) 


Compressibility Equations 

There are two volumes in the fuel-control valve: (1) the volume between the pressure 
reducing orifice and the control orifice and (2) the volume between the control orifice and 
the valve output. 

If the fuel is considered to be a compressible fluid, the dynamic relation between 
pressure and volumetric flow through a control volume is as follows (ref. 2): 


p _ ^ vol, in ^vol, out 
0 " V 


( 8 ) 


For the control volume between the reducing valve orifice and the control orifice: 


^vol, in ^vol, out ^r ” ^c + A r x r 


(9) 


and 


V=V r“ A r x r (10) 

The third term on the right side of equation (9) is the fuel flow generated by the reducing 
valve piston motion. Substituting equations (9) and (10) into equation (8) yields 


P r _ q r - q c + A r x r 
0 ' ‘v r - Vr 


(ID 


If, for all x r , V r » A r x r , then equation (11) reduces to 

b q„ - q„ + A„x„ 

*r _ ^r c r r 



For the control volume between the control orifice and the output 


( 12 ) 
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(13) 


^vol, in ” ^vol, out 


= % - % ~ Vr 


and 


v = v c 


+ Ax 
r r 


Substituting equations (13) and (14) into equation (8) yields 


Pc 

P 


Vr 


(14) 


(15) 


for V c » A r x r . 


Force-Balance Equations 

Fuel pressure actuates the reducing valve piston. Therefore, a force-balance equa- 
tion describes the motion of this spring mass system. For the valve being evaluated in 
this report, the area of both faces of the reducing valve piston are equal. The equation is 

M r* r + B r x r + k r x r = A r^ p c " P r ) + F (16) 

where B r is a friction force or a force due to velocity. Its value could not be measured 
and was chosen to give a damping of 0. 2 for the left half of equation (16). 


Lossless Transmission Line 

In a typical test installation, the fuel-control valve may be located some distance 
from the engine load. If this is the case, then some type of fuel line will connect the 
valve to this spray nozzle load. For the purposes of this report, this is assumed to be a 
lossless hydraulic transmission line. 

In reference 3 two equations that relate the flows and pressures for the line are of 
use in this study. These are, in the time domain 

P L (t) = p c (t - o) + Z fq L (t) + q Q (t - cr)j (17) 


6 



q o (t) = q L (t - cr) +1 [p c (t) - p L (t - a)] 


(18) 


Nonlinear Block Diagram 

The nonlinear block diagram (fig. 3, p. 11) represents symbolically the function of 
equations (2), (4), (6), (12), and (15) to (18). The inputs to the block diagram are p^, 
p g , F, and x c< 


LINEARIZED FUEL VALVE REPRESENTATION 

In actual usage for a dynamics research program, the fuel-control valve is used to 
modulate fuel flow dynamically about some steady-state operating level of fuel flow. 

Thus, valve operation will involve small deviations about a selected operating point. This 
type of operation makes linearized analysis of the describing nonlinear equations an ideal 
method for evaluating the valve small-disturbance performance. Appendix B includes the 
linearization of the basic valve equations. Since most experimental modulating signals 
will be sinusoidal in nature, the linearized equations have been converted to the Laplace 
domain. Also, from experimental observation and from the nonlinear simulation, the 
pressure-reducing valve closely regulates the control orifice pressure difference 
(P r - P c ) to frequencies much higher than those of interest herein. This feature simplfies 
the resulting linearized equations. The equations resulting from these manipulations and 
the preceding observations are stated as follows: 


where 


AQ C = K x AX C 


r r 


A r (V r + V c ) 


K 3 £ 


sK(^ + ^ 


j8 K c 


S > AP C = AQ c - U + -L SjAQ c 


AQ t = K- AP, 


K 

K 


1 

3 




t=T 


t=T 


(19) 

( 20 ) 

( 21 ) 

( 22 ) 

(23) 
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(24) 


k 4 = 


K x 
r r 


2 Vp S -Pi 


lt=T 


K t 


K 5 = 




t=T 


(25) 


Equation (20) defines the output characteristic of the fuel-control valve. This equa- 
tion can be used to determine the transfer functions of this device when operating with any 

load whose input flow AQ is linearly related to the valve output pressure AP . 

o ^ 


Closed Coupled Load 

If the load orifice is connected directly to the fuel- control valve outlet, then 
AP^ = AP c and AQ q = AQ^. Substituting equation (21) into (20), the following transfer 
function is obtained: 


AQl 

AQ C 


[ VV * Vc V +\L+ : 

L K 3 K 6 (! J Lk 5 /3 K 3 K 5 J 


(26) 


Equation (26) is in the general form of a second-order transfer function with the natural 
frequency and damping defined as follows: 

(1) For the natural frequency 




K 3 K 5 p 


A r (V r + V c ) 


(27) 


(2) For the damping coefficient 


lK + A r< K 4 + K 5? 1 ^ 

2b K 3 J \K 5 A r (V r + V c ) 


(28) 
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Transmission Line Coupled Load 


Equations (17) and (18), which describe the long fuel transmission line, can be 
linearized directly and result in 


Ap L (t) = Ap c (t - a) + Z [-Aq L (t) + Aq 0 (t - a)J 

(29) 

Aq 0 (t) = Aq L (t - a) + — [Ap c (t) - Ap L (t - a)] 

z 

(30) 


Application of the Laplace transform converts equations (29) and (30) into 

AP l = AP c e" S(T + Z (-AQ l + AQ o e' Scr j (31) 

c " AP L e ' Sa ) < 32 > 


AQ q = AQ L e _Sa + 


i( AP 


When the long fuel transmission line is used to connect the valve to the load orifice, 
AP c is no longer equal to AP L and AQ q is no longer equal to AQ l> Therefore, equa- 
tions (20), (21), (31), and (32) can be combined by standard algebraic means and result in 
the following transfer function: 



(33) 


If s = jco, equation (33) becomes 


AQ X 

AQ„ 


cos uicr + 


1 ^ A r< V r + V c> Z 
ZK r K„/3 


sin co< 



A r< V r + V c> 


k 5 k 3 ^ 


(34) 


cos too’ - sin o>a| 


where 


K, 


K 

1/3 


A r< K 4 

Ko 


+ K 5>] 1 

. Kc 


(35) 


9 



and 


K 


7 “ 



(36) 


-1 J 
When Z = Kg , the load impedance is matched to the transmission line impedance, and 

the amplitude frequency response of the fuel system will show no degradation due to line 
resonances. However, matching the line impedance to the load will not reduce the dead 
time caused by the transmission line. A FORTRAN IV program for computing the fre- 
quency response of equation (34) appears in appendix C. 


COMPARISON OF DYNAMIC PERFORMANCE, LOAD CLOSE COUPLED 
Simulated Nonlinear Model Performance 

The nonlinear block diagram of figure 3 was simulated on an analog computer using 
the parameter values of table I. Figure 4 shows the fuel-flow frequency response for 
this model with a close coupled load to a small control orifice area input signal. The 
magnitude of this input signal corresponds to 5 percent, peak-to-peak, of maximum fuel 

O 

flow, where maximum fuel flow is 38. 5 cubic inches per second (630. 9 cm /sec). 


Experimental Performance 

Figure 5 is a photograph of the experimental valve, disassembled to show the ar- 
rangement of its internal parts. The actual frequency response data for this configura- 
tion is shown in figure 6. The nonlinear analytical results of figure 4 have been over- 
layed for comparison purposes. There is close correlation for all values of fuel flow. 
Experimental fuel flow was determined by measuring the pressure drop across the load 
orifice. 


Linearized Model Performance 


For an input of 5 percent, peak-to-peak, of maximum fuel flow, 
sponses of the linear and nonlinear models are found to be identical. 


the frequency re- 
Therefore, the 
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Figure 4. - Fuel-flow frequency response from simulated nonlinear model. Close 
coupled load. 


TABLE I. - FUEL CONTROL VALVE MECHANICAL CONSTANTS 


2 2 

Reducing valve piston area, in. ; cm 

1.0; 6.45 

Reducing valve friction, B, (lbf)(sec)/in. ; (N)(sec)/cm 

0.0633; 0.1109 

Reducing valve spring bias, F, Ibf; N 

140; 623 

Control orifice constant, K c , in. V(sec)(lbf 1//2 ); cm^/(sec)(N 1 ^) 

16.25; 126.26 

Load orifice constant, K L , in. 4 /(sec)(lbf* ,//2 ); cm 4 /(sec)(N 1//2 ) 

1.397; 27.57 

Reducing valve orifice constant, K r , in. 3 /{sec)(lbf 1//2 ); cm 3 /(sec)(N 1//2 ) 

70.4; 547.0 

Reducing valve spring constant, k r , Ibf/in. ; N/cm 

100; 175 

2 2 

Reducing valve piston mass, M r , (lbf)(sec )/in. ; (N)(sec )/cm 

2.51X10" 4 ; 4.40 

2 

Fuel supply pressure, p g , psi; N/cm 

650; 448 

3 3 

Control orifice volume, V , in. : cm 

c ? 

7.5; 122.9 

3 3 

Reducing valve volume, V , in. ; cm 

5.55; 90.95 

5 5 

Hydraulic line surge impedance, Z, (lbf)(sec)/in. ; (N)(sec)/cm 

9.235; 16.18 

2 

Fuel bulk modiolus, /3, psi; N/cm 

1.5X10 5 ; 1.04X10 5 

2 4 2 4 

Fuel mass density, p, (lbf)(sec )/in. ; (N)(sec )/cm 

7. 5X10" 5 ; 8.02X10" 6 

Hydraulic time delay, a, sec 

2. 6x10" 3 

2 2 

Hydraulic line area, a, in. ; cm 

0.36317; 2.343 

Hydraulic line length, l , in.; cm 

120; 305 




Figure 6. - Comparison of experimental and analytical fuel-flow frequency response. 
Close coupled load. 


linear model is quite adequate for dynamic performance evaluations at these relatively 
low levels. It should be noted that the valve natural frequency and damping (eqs. (27) 
and (28)) are dependent on several fixed valve parameters as well as several lineariza- 
tion constants. The linearization constants vary with the steady-state operating point. 
This variation is evident in the frequency response curves of figures 4 and 6 and in equa- 
tion (26) for different levels of steady- state fuel flow. It should be noted that dynamic 
performance or bandwidth degrade as the steady-state fuel flow operating point increases. 
These three steady- state conditions are tabulated in table n. 
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TABLE D. - FUEL CONTROL VALVE OPERATING CONDITIONS 


q o 

[Maximum fuel flow, 38. 5 in. /sec (10. 0 gal/min or 630. 9 cm /sec); 
fuel supply pressure, 650 psi (448 N/cm 2 ). ] 


Percent 

maximum 

flow 

Reducing 

valve 

pressure, 

Pr 

Control 

orifice 

pressure, 

Pc 

Control orifice 
length, 

x c 

Reducing valve 
length, 

x r 


psi 

N/cm 2 

psi 

N/cm 2 

in. 

mm 

in. 

mm 

31 

213 

147 

73 

50 

63x10* 3 

1.57 

8. 13x10* 3 

0.206 

54 

362 

250 

222 

153 

109 

2.77 

17.38 

.441 

71 

523 

361 

383 

264 

144 

3.66 

34. 15 

.867 


DYNAMIC PERFORMANCE WITH A LONG LINE 
Performance of Simulated Model 


The nonlinear block diagram of figure 3 was simulated on an analog computer using 
the parametric values of table I. The 2. 6-millisecond time delay was simulated with a 
fourth-order Pade approximation (ref. 4). Dynamic performance was evaluated for the 
case of a 10-foot (3.048-m) long, 0. 75-inch (1.905-cm) outside diameter, stainless- steel 
fuel transmission line at 54-percent maximum steady-state fuel flow using the same load 
orifice as the close coupled case. The performance for this condition is given in figure 7. 
The performance of the close coupled nonlinear model at 54-percent maximum fuel flow 
is overlayed for demonstration purposes. The effect of the transmission line is consid- 
erable, degrading the system bandwidth by nearly a 4 to 1 margin. 



Figure 7. - Comparison of nonlinear model fuel -flow frequency responses for close 
coupled and transmission line coupled load. Steady-state fuel flow , 54 percent 
maximum. 
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Experimental Performance 


The experimental frequency response curve obtained at 54-percent maximum flow 
condition is shown in the curve of figure 8. The analytical performance of figure 7 as ob- 
tained from the simulated nonlinear model, has been overlayed for comparison purposes. 
These results were obtained with a 10-foot (3. 048-m) long, 0. 75-inch (1. 905-cm) outside 
diameter line. The theoretical effect of line resistance in this particular case was 
checked (ref. 5) and found to make the model results fall closer to the experimental data. 
However, the improvement was not significant enough to require the modification of the 
basic flow equations (17) and (18). In installations in which smaller diameters or longer 


t/1 o 

CL) 

0 



O 



Figure 8. - Comparison of experimental and analytical fuel -flow frequency response. 
Steady-state fuel flow, 54 percent maximum; transmission line coupled load. 


lengths of transmission are used, the effect of this theoretical line resistance may not be 
negligible. 

It should be noted that the attenuation of the analytical model performance is not as 
great as that of the actual experimental performance. The curves do, however, possess 
the same ?t valley" frequency or frequency at which attenuation is a maximum. No ex- 
planation for the discrepancy is available. The model, however, should be adequate to 
approximately predict the degradation that will occur at various frequencies due to the 
addition of the long transmission line. 


Long Line Linear Model 

Figure 9 gives the frequency responses for the three conditions of fuel flow as tabu- 
lated in table II. The performance of the linear model is again identical to the nonlinear 
model performance. 
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Frequency, Hz 


Figure 9. - Frequency response of linear model at different fuel flows. Transmission 
line coupled load. 

Figure 9 shows how the fuel-flow frequency response varies at the three levels of 
steady-state fuel flow. As mentioned earlier, when the line impedance is equal to Kg* 
the fuel-flow amplitude frequency response shows no degradation due to the transmission 
line. When the product ZKg is less than 1, a valley will occur in the fuel-flow ampli- 
tude frequency response at a frequency somewhat less than l/4or hertz, and, when ZKg 
is greater than 1, a peak will occur. Thus, the product ZKg could be used as a meas- 
ure of the line to load matching 

(37) 

a 2 Vf^~ 2aq L 

For a given steady-state flow condition, equation (37) shows that the impedance match is 
proportional to the square of the load orifice constant (thus, the square of the orifice 
area) and inversely proportional to the transmission line cross-sectional area. Equa- 
tion (37) also shows that, if the steady-state fuel flow is changed, the impedance match is 
changed. Thus, if the line impedance is matched to the load for one level of fuel flow, 
the line and load will not be matched at other levels of fuel flow. For the experimental 
system described in table I, equation (37) is less than one for all three levels of fuel flow. 

DISCUSSION OF DYNAMIC PERFORMANCE AND RECOMMENDATIONS 

From the response curves of figure 5, it may be concluded that the nonlinear analyt- 
ical model is a very good representation of the fuel system with close coupled load. The 
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second-order linear transfer function describing this configuration also accurately repre- 
sents the dynamic system. The main advantage of the transfer function linearization is 
that the fuel system flow response for various conditions can be predicted almost by in- 
spection. 

For example, different supply pressures and/or different load orifices may be used. 
Both parameters affect the values of u> n and 6, the natural frequency and damping of 
the valve. Figures 10 and 11 show how ci> n and 6 vary, respectively, as functions of 
the supply pressure with a specific load orifice constant. The value of is 1. 397 
inch^ per second per pound"* - ^ (27. 5 cm V(sec)(N^^)). We can see that, as the supply 
pressure is decreased, o> n decreases and damping increases. Both effects will degrade 
the dynamic performance capabilities of the fuel valve. Thus, to operate unattenuated 
to the highest frequencies, supply pressure should be maintained at the higher level 
(650 to 800 psi (450 to 550 N/cm 2 )). 

Figures 12 and 13 show how u> n and 6 vary, respectively, as functions of the load 
orifice constant K T with supply pressure held constant. The value of p„ is 650 psi 
(448 N/cm^). The natural frequency decreases as the load orifice constant decreases 
(load orifice area decreasing). However, for the three flows considered, the damping 
remains nearly constant (between 0. 8 and 1. 2) over a fairly wide range of load orifices 
(1.65 to 3.8 in.^/(sec)(lb'^) or 32 to 74 cmV(sec)(N^ 2 )). 

The long line linear model is also considered to be sufficiently accurate for determin- 
ing fuel system flow response with a lossless transmission line coupling the load. The 
transmission line degrades the amplitude response of the fuel system flow (fig. 7) as well 



Figure 12. - Effect of load orifice constant on valve natural frequency. Supply 
pressure constant at 650 psi(448 N/cm 2 ). 
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Load orifice constant, K|_, cm^/isecKN 1 2 ) 


Figure 13. - Effect of load orifice constant on valve damping. Supply pressure 
constant at 650 psi (448 N/cm 2 ). 

as adding phase shift due to its time delay. It is recommended that the use of a transmis- 
sion line be avoided so that the user can take full advantage of the flow response capabil- 
ities of the fuel-control valve. If the transmission line must be used, the analytical mod- 
els developed in this report can be used to approximately predict the effect of the trans- 
mission line on the system flow response. 


CONCLUDING REMARKS 

An analysis of the fuel-control valve with a pressure reducing regulator, including 
the effects of a long hydraulic transmission line has been presented. Nonlinear and lin- 
ear dynamic models of the valve were developed for both the close coupled load and the 
load coupled by a long hydraulic transmission line. These models predict the perform- 
ance of these systems up to 100 hertz. Therefore, these models can be used successfully 
to predict fuel-flow dynamic performance for different configurations and various steady- 
state operating points. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, April 26, 1969, 

720-03-00-69-22. 
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APPENDIX A 


SYMBOLS 

..22 

A area, in. ; cm 

2 2 

a hydraulic transmission line cross-sectional area, in. ; cm 
B sliding friction, lb-sec/in.; (N)(sec)/cm 

C orifice discharge coefficient 

c velocity of sound in A-l jet fuel, VlVp, in. /sec; cm/sec 
F spring bias force, lb; N 

h orifice width, in. ; cm 

i \PI,j 2 = -i 

K with alphabetic subscripts, orifice constant; with numeric subscripts, lineariza- 
tion constant 

k reducing valve spring constant, lb/in.; N/cm 
l length of hydraulic transmission line, ip. ; cm 

9 9 

M mass of piston, (lb)(sec )/in. ; (N)(sec )/cm 

9 

AP Laplace transform of Ap, psi; N/cm^ 

9 

p gage pressure, time domain, psi; N/cm 

o o 

AQ Laplace transform of Aq, in. /sec; cm /sec 

o o 

q volumetric flow, time domain, in. /sec; cm /sec 
S Laplace transform operator, sec"* 

T time at which steady-state values exist 
t time, sec 

3 3 

V volume downstream of variable orifice, in. ; cm 
AX Laplace transform of Ax, in. ; cm 
x variable orifice length, time domain, in. ; cm 

Z hydraulic transmission line surge impedance, Vpfi/ a, (lb)(sec)/in. (N)(sec)/cm^ 

p 

/3 bulk modulus of A-l jet fuel, psi; N/cm 

5 second- order transfer function damping coefficient 

d phase angle, deg 
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p 


2 4 2 4 

mass density of A-l jet fuel, (lb)(sec )/in. ; (N)(sec )/cm 
a hydraulic transmission line time delay, l / c, sec 

it) frequency, rad/ sec 

Subscripts: 
c control 

d drain 

L load 

n second order transfer function natural frequency 

o fuel valve outlet 

r reducing valve 

s supply 

vol, in volume in 

vol, out volume out 

1-7 linearization 
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APPENDIX B 


LINEARIZATION OF FUEL VALVE EQUATIONS 

Repeated are the nonlinear differential equations for the case of the fuel control 
valve and load orifice. 


q c = K c x cVP r -P c 

(4) 

q r = K r x r Vp s - p r 

(6) 

q L = K L 

(2a) 

P r _ q r - q c + Vr 
(3 V r 

(ID 

Pc _ «■: - % - Vr 

f* V c 

(15) 

M„x„ + Bx + k x = A (p„ - p ) + F 
r r r r r r r vr c r r ' 

(16) 


From experimental observation and from the nonlinear simulation, the pressure re- 
ducing valve closely regulates p r - p c to frequencies much higher than those of interest 
here. Thus, equation (16) becomes 

0 = A r (p c - P r ) + F (Bl) 

Applying the Laplace transform, equation (Bl) becomes 

AP c = AP r (B2) 

Equation (4) is linearized as follows: 

Aq c = K : Ax c + K 2 Ap r - K 2 Ap c (B3) 


where 
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( 22 ) 


K 1 = K c VP r ' P c 


I t=T 


and 


K 2 = ~ K c x c / 1 = 

2 VP r " Pc 


(B4) 


t=T 


Applying the Laplace transform and substituting equation (B2) yields for equation (B3) 

AQ c = K x AX C (19) 

In the same manner, equations (6) and (2a) are linearized as follows: 

Aq r = Kg Ax r - K 4 Ap r (B5) 


where 


K 3 = K r V 'Ps " p r 


t=T 


(23) 


and 


K 4 = i K„x 


2 ^Vi^r 


1 t=T 


Aq L = K 5 Ap L 


(24) 


(B6) 


where 


K = i K -J— 
5 2 L 


t=T 


Applying the Lap lac e transform yields for equations (B5) and (B6) 

AQ r = Kg AX r - K 4 AP c 


(25) 


(B7) 


where 
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AP r = AP c (B8) 

and 

aq l = k 5 AP l (21) 


Similarly, equations (11) and (15) yield 


where 


AP c S-i-= AQ r - AQ c + SA r AX r 


AP = AP 
r c 


(B9) 


AP c S-^=AQ c -AQ 0 -SA r AX f 


P 


(BIO) 


By standard algebraic methods, equations (B7), (B9), and (BIO) can be reduced to 
yield the following equation: 


C- 


•A r (v r + V c r 


s 2 + 


L k 3 ^ 




( 20 ) 


Equation (20) defines the output characteristic of the fuel-control valve. This equation 

can be used to determine the transfer function of this device when operating with any load 

whose input flow ACD is linearly related to the valve outlet pressure AP . 

o ^ 
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APPENDIX C 


DIGITAL COMPUTER PROGRAM AND FLOW DIAGRAM 


DIMENSION F< 25) ,W(25) ,C(25> , S(25 ) 

REAQ{ 5,2) P$,Z,TE, AKL , AKR , VC , VR , B, AR 

WR I TE ( 6,2) P S,Z, TE , AKL , AKR, VC, VR ,6 ,AR 

D3 I [-1,25 

EF 1=1- 1 

EF 2 = EF 1/ LO.O 

F ( I )=10.0**EF2 

W ( I ) -6 .2d3LB5*F ( I ) 

C< I ) = CGS{ W ( I )*TE) 

1 S( I ) = S INK W(I ) * TE ) 

3 REALM 5,2) PR ,PL , XR 
SQ = SQR 1 ( P S-P R ) 

AK3=AKR*SQ 

AK4-0 . 5*AKR* XR/SQ 
A < 5 = 0 • 5*AKL/SQRT(PC) 

A< 6 — ( VC/3 +AR* ( AK4+AK5) /AR3) / AK5 

A< 7 = { VC/S +AR * ( AK4+1.0/ (AK5*Z*Z) ) /AK3 ) *Z 

CK 1 = L • 0/{ Z*AK5) 

C< 2=AR*( VR + VC) *Z /( AK3*B) 

CK 3=A R*( VR+VC) /(AK3*AK5*3) 

WR I TE ( 6,4) 

OJ 5 1-1,25 

TI=A<6*W( I)*C(IH-CCia-:K2*hl(I)*hi([))*Sm 
TR — AK 7*W ( I )*$< I ) + ( 1. 0-CK3*W(I ) *W( I))*C( I ) 

GS = T I # TI + TR^ TR 

IF (GS .LT. l.OE- 10) SO TO 10 

AM AG-SQ-R T ( 1. O/GS) 

I F ( TR ) 7,8, 9 

7 PHASE — LBO.O +5 7.2 9 577*ATAN(-TI/TR) 

GJ TO 6 

8 I F ( T I ) 11,11 ,12 

1 1 PHASE — 2/0.0 

GJ TJ 6 

12 P HA SE =— 90 . 0 

GO TO 6 

9 PHASE- 57. 2 95 77*A TAN( - TI/TR) 

6 WR ITE ( 6, 1 3 )F (I) , AMA3 , PHASE 

GJ TJ 5 

10 WR ITE (6,14) 

5 CJHTIN1UE 
GJ TJ 3 

2 FORMAT ( 3E 2 0. 8) 

4 FORMAT ( 1H 1 , 1 OX , 9HF RE QUENC Y , 1 LX, 9HMAGNI TUDE ,1 LX ,5H PHASE//) 

13 FORMAT ( IX, IP 3E20. 8) 

14 FJRMAT (IX, 12HMAG TO LARGE) 

EMO 




o 
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